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, NOTATIONS 

A — the right ascension of the satellite, 

r — the distance of the satellite from the center of the Earth, 
r° — the unit vector directed from the center of the Earth toward the 
satellite, 

k, i u,v — the equatorial components of r° , 
g — the mean anomaly of the satellite, 
f — the true anomaly of the satellite, 
co — the argument of perigee of the satellite, 

— the right ascension of the ascending node of the satellite, 

77 = w+ A — the longitude of the perigee of the satellite, 

L = g + 77 — the mean longitude of the satellite, 

i - the inclination of the orbital plane of the satellite relative to the 
equatorial plane of a given epoch, 
e — the eccentricity of the satellite's orbit, 
a — the semi-major axis of the orbit of the satellite, 
n — the mean motion of the satellite, 

R — the unit vector normal to the orbital plane of the satellite, 

SL, 877, SA,Si - the tidal perturbations in L, 77, A and i, respectively, 
m' — the mass of the Moon (Sun), 

A' — the right ascension of the moon (Sun), 
x*, y', z' — the equatorial rectangular coordinates of the Moon (Sun) referred 
to the mean equator and equinox of a given date, 
r' — the distance of the Moon (Sun) from the center of the Earth, 
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A.' , fj.' , v' — the equatorial components of the unit vector r'° directed 
from the center of the Earth toward the Moon (Sun), 
a* - the semi-major axis of the lunar (solar) orbit defined in such a 
manner that the constant part in the expansion of a'/r 1 be equal to one, 
G — the gravitational constant, 

M - the mass of the Earth, 

R — the equatorial radius of the Earth, 
e — the eccentricity of the meridian of the Earth, 
k , k , • • • — Love numbers, 
a = R/a, 

a' = R/a* - the lunar (solar) parallactic factor, 
ft — the tidal disturbing function, 


— — ■ — the component of the disturbing force normal to orbital plane 
o Z 

of the satellite. 


P nm — associated Legendre functions, 
I - the idemfactor. 
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INTRODUCTION 

In the present article we develop and extend the ideas on the tidal perturba- 
tions of artificial satellites formulated in a previous article by Musen and 
Estes (1971). The magnitude of the tidal effects depends upon the elastic properties 
of the Earth as described by Love numbers. The Love numbers appear as the 
coefficients in the expansion of the exterior tidal potential in terms of spherical 
harmonics (in Maxwellian form). In the previous article, we considered the 
Love numbers as constants over the whole Earth and substituted their "global" 
values into the potential and into the differential equations for the tidal perturba- 
tions in the motion of satellites. 

These global values can be understood to be doubly averaged local Love 
numbers, the averaging process being performed over the whole surface of 
the Earth. 

Such an approach yields many interesting and important results (Smith, 
Kolenkiewicz and Dunn, 1971). However, there are discrepancies between the 
values of Love numbers as obtained from the observations of satellites. This 
suggests a necessity for further sharpening the idea and making use of a single 
averaging process performed only along the parallels of latitude. In this respect, 
we follow in the footsteps of Kaula (1969). 

This single averaging process preserves the additional long period tidal 
effects (with periods of a few days or more) which were omitted in the previous 
work. It also eliminates the short period effects with periods of one day or less. 

It seems that only the k 2 and k 3 Love numbers and their deviations from con- 
stants are of any observational importance. The effects associated with k 4 . 
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k 5 , . . . will be damped out by the presence of higher powers of the lunar 
(solar) parallactic factor. 

It must be emphasized that if the deviations of k 2 and k 3 from pure constants 
on the surface of the earth are not taken into consideration, then the danger arises 
that the perturbations which are caused by the variability of k 2 and k 3 will be 
misinterpreted as the perturbations caused by k 4 , . . . Then, the wrong values 
of k 4 , . . . would be obtained. 

In the present article, we obtain the expansion of the exterior tidal potential 
using the theorem of Dirichlet. We also obtain.the differential equations for 
tidal perturbations of elliptic elements taking into account the deviation of k 2 
and k 3 from constants along the parallels of latitude. 

These differential equations are in a form convenient for numerical inte- 
gration as well as for the trigonometrical solution (Musen and Estes, 1971), with 
the arguments " , fl of the satellite and with the standard arguments M\f, 

D, and r of the lunar theory. The short period terms with periods equal to 
the period of revolution of the satellite or less, were removed from the differ- 
ential equations by averaging over the instantaneous orbit of the satellite. 

A numerical integration and the prediction of the tidal effects is feasible 
for an interval of a few months. 

Using this method, a general idea about very long period influences can be 
obtained. For example, the influence caused by the rotation of the lunar node 
can be detected by integration over a number of years. 

A full collection of formulas is given to facilitate the programming. 
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DISTURBING FUNCTION 


We make use in the present article, similarly as we did in the previous 
work (Musen and Estes, 1971), of the Maxwellian form of the expansion of the 
tidal potential into a series of products of spherical harmonics. The first 
factor in each product is a polynomial in the components V , // , v' of the 
unit vector directed from the center of the Earth toward the Moon, and the second 
factor is a polynomial in the analogous quantities k, j±, v for the satellite. 

We resort to this representation of the disturbing function, and not to the 
representation in terms of the elements of the motion of the Moon, because at 
the present time we do not possess the expansion of the lunar osculating elements 
into a periodic series with the standard arguments -t, V , F and D of the lunar 
theory. The raised standards of the accuracy of computations do not favor the 
substitution of the lunar mean elements into the differential equations for the 
lunar or tidal perturbations of the satellite. 

We know, however, the expansions of the rectangular coordinates of the 
Moon, and consequently also of k' , /j.' , v' , and their numerical values with an 
accuracy which permits us to determine the tidal effects of even higher order. 

The tidal disturbing function on the surface of the Earth has the form 




A' a. 


i j 





( 1 ) 


where we set 



Ai = t i-(l-3v' 2 ) = -ip 20 (O. 

a; = +|(\' 2 -m' 2 ) = + |p 22 O'') 20524 '' 

a; = t 3 \ V =+iP 22 (v')sin2A\ 

A' = + 3 v* = P 21 (*') cos A', 

A; = + 3 ft v‘ = P 21 (v') sin A'. 

B! = +1\' (1 - 5 v' 2 ) = --jP 31 (v*) cos A', 

B' =+|/ (1 -5 v' 2 ) = -iP 31 (v')sinA', 

b; = +ii/' (3-5 v' 2 ) = -ip 30 O''), 

b; = +|A* (A' 2 -3/ 2 ) = +^ P 33 O'*) cos 3A '’ 

b; = + |/ (3 X.' 2 -m' 2 ) = + ^4 P 33 OO sin3A '’ 
B ; = + ^v' (V 2 -/x' 2 ) = + ip 32 O'') cos 2 A', 
B' 7 = + 15 K' it v' = +|p 32 O'') sin 2A '- 
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a 0 = k 2 r 2 (1 - 3 v 2 ) = - 2 k 2 r 2 P 20 (v); 


a 1 = k 2 r 2 (A . 2 - /j?) = + | k 2 r 2 P 22 (v) cos 2 A, 

a 2 = k 2 r 2 = + i k 2 r 2 P 22 (v) sin 2A, 

a 3 = k 2 r 2 X y = + i k 2 r 2 P 21 (v) cos A, 

a 4 = k 2 r 2 /iv = + ^ k 2 r 2 P 21 (v) sin A. (3) 

b t = k 3 r 3 \ (1 - 5 v 2 ) = - k 3 r 3 P 31 (v) cos A, 

o 

b 2 = k 3 r 3 /Lt (1 - 5 V 2 ) = - k 3 r 3 P 31 (v) s in A, 

b 3 = k 3 r 3 v (3 - 5 v 2 ) = - 2 k 3 r 3 P 30 ( v ) , 

b 4 = k 3 r 3 \ (\ 2 " 3 M 2 ) = + -jj k 3 r 3 P 33 (v) cos 3 A , 

b 5 = k 3 r 3 ft (3 X 2 - fj?) = + -L k 3 r 3 P 33 (v) sin 3 A, 

b 6 = k 3 r 3 v (\ 2 ~ M 2 ) = + j^k 3 r 3 P 32 (v) cos 2 A, 


b ? = k 3 r 3 Kfiv= + k 3 r 3 P 32 (v) sin 2 A. 


( 3 *) 
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and k 2 , k 3 , . . . are the Love numbers. They characterize the local elastic 
properties of the Earth on its surface. The Love numbers are usually assumed 
to be constants and their so-called "global" values are substituted into the ex- 
pansion of the tidal disturbing function. The "global" Love numbers represent 
the averaged elastic properties of the Earth as they are being observed by the 
satellite from outer space. The introduction of the global Love numbers yields 
interesting and important results confirmed by observations. All existing 
analytical expansions or the numerical integration of the tidal effects consider 
only k 2 . 

The inclusion of the k 3 -part, and possibly the higher order parts, represents 
a more subtle problem and attention must be given to some local properties of 
Love numbers and also to the deviation of the Earth from a sphere. 

We assume, together with Kaula (1969), that the Love numbers are functions 
smoothly varying on the surface of the Earth and can be expanded into a series 
of spherical harmonics. We do not know yet the values of coefficients in this 
expansion, but we can expect that in our problem, because of the averaging 
process performed by the rotation of the Earth, the dependence of k 2 , k 3 , ... on 
the latitude will be stronger than dependence upon the longitude and we can set 
(Kaula, 1969): 

k l =k U + hl P 10 (*0 +k ^ 2 P 20 O') + ■ 

The neglected terms in this expansion are periodic with a period of one day or 
less. At the present time they can be considered as short periodic and therefore 
can be omitted. 
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Thus, following Kaula (1969), the averaging process will be performed only 
along the parallels, and not over the whole Earth. The values of k 2 as obtained 
from the observation of satellites by different authors (Kozai, 1965), (Newton, 
1968), (Smith, et. al., 1971), (Musen and Estes, 1971) are not in complete agree- 
ment. The determination of k 3 from the motion of satellites still must be done. 

In order to facilitate the selection of the numerically significant part of the 
tidal potential we computed estimates of the factors attached to the terms in the 
expansion of this potential. We have the following information on the basis of 
present knowledge: 


k 2Q ~ 0.3, k 3Q ~ 0.1, k 22 ~ 0.1 

and we believe that we are on the safe side if we assume the following preliminary 
estimates: 


k 21 ~ 0.01, k 31 ~ 0.01, k 4 ~ 0.1 

Using these numbers we obtain the following estimates for the indicated factors: 


1°. direct k 2 part: 


2°. direct k 3 part: 


3°. direct k 4 part: 


m' 

¥ 


a' 3 k 20 ~ 1.7 x 10~ 8 , 


m' 

¥ 


a' 4 k 30 ~ 1 x IQ -10 


— a' 5 k. ~ 1.6 x 10“ 12 
M 4 


7 



4°. terms representing the influence on 2° of the deviation of k 2 from a 
constant: 

— a' 3 k ~ 5.6 x 10" 10 

M 21 

5°. terms representing the influence on 3° of the deviation of k 2 from a 
constant: 

— a' 3 k„ - 5.6 x 10- 9 , 

M 22 

6°. terms representing the influence on 3° of the deviation of k 3 from a 
constant: 

— a' 3 k „ - 5.6 x lO' 10 

M 31 

7°. terms representing the influence of the oblateness on 2°: 

ilL a' 4 k on e 2 - 1.9 x 10* 12 

M 20 

8°. terms representing the influence of the oblateness on 3°: 

— a' 3 k, ft e 2 - 1.2 x lO" 10 

M 20 

We see from this table that the part 3° associated with k 4 constitutes only 1% 

2 

of the terms 8° having the same trigonometrical arguments as 3° and k 20 e as 
a factor. 
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Also, we maintain that the terms having k s (s > 3) as a factor will be less 
influential than the terms having the same trigonometrical arguments and factors 
of the type k 20 x k 2j , k 30 x k 3j . 

In this case it will be difficult to determine k s (s > 3) Love numbers. The 
parts of the tidal potential having k 2 and k 3 as factors will be the parts which 
have dominant physical and numerical significance in the theory of tidal 
perturbations of close satellites. ■' 

For this reason, we include in the disturbing function terms of the types 
1°, 2°, 4°, 5°, 6° and 8°, and omit the terms 3° and 7°. 

Substituting 

k 2 ~ k 20 + k 21 P 10 "*■ k 22 P 2 0 
k 3 =k 30 +k 31 P 10< V ) + 

into (3) and (3 ’) , retaining only the significant terms and redeveloping the results 
in terms of spherical harmonics in form convenient for the continuation of the 
potential into outer space, we obtain 

a ° = r2 ["I 1 * 22 ‘I 1 " 21 V + ( k2 ° + I k22 ) (1 ” 3v2) 

+ -| k 2i V (3 - 5 V 2 ) k 22 (3 - 30 V 2 + 35 v 4 )J , 
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I 1 =r 2 (\ 2 -M 2 ) £ + ( k 20 - 7 k 2 2 ') +k 


21 * - 4 k 22 C 1 - 7 ^ 2 ) 


i- 


a - r 2 


a 3 = r* 


a 4 = r- 


[+ ( k 20 - 7 k 22 ) + k 21 ^ ^ k 22 C 1 - 7 ^ 2 )] ■ 
X \_ + J k21 + ( k20 + y k 22 ) V -^ k 21 ^ - 5y2 > 

H4 

M [+| k 21 + ( k 20 + y k 2 2 ) v - | k 21 C 1 - 5v2 > 


-14 k22 " 


M4 


and 


( 4 ) 


b i r ' = r 3 ^ y k 31 ^'+ k 30 (1 - 5 V/ 2 ) + y k 31 V |l - I v 2 jj , 

b 2 - 1-3 M [-y k 31 V +k 30 ( l - 5^ 2 ) + yk 31 v (l 

b 3 ■ = r 3 |^+| k 3 x (1 - 3v 2 ) +k 3Q v (3 - 5v 2 ) -ik 31 (3- 30v 2 + 35v 4 )J, 

b 4 = r 3 k (X. 2 - 3 fj?) (k 30 + k 31 v ), 

b 5 = r 3 /x (3 X. 2 - fj. 2 ) (k 3Q + k 31 v), 

b 6 = r 3 (4 2 - fj?) |\ Ik 31 + k 30 1 / - Ik 31 (1 - 7 i^ 2 )J , 

b 7 = r 3 \/x j^+ j k 31 + k 30 v - I (1 - 7 y 2 )J • ( 4? ) 
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For the extensions of a . and b . into outer space we have with the same 
accuracy as before: 





a' 3 1 3 

a a J 

a 



v 


+ 





(1-3 z, 2 ) 



,'3 



(3 v - 5 v 3 ) 





a 


a 


'3 



(3 - 30 i^ 2 + 35 v 4 ) , 





a 


a 


'3 




+ k 


,'3 


21 


a 3 a' 3 



(A 2 - yu 2 ) 




£la 4 a' 3 (l) 5 (X 2 -^ 2 ) (1-7 i/ 2 ), 



+ k,.^a 3 a' 3 
a 



+ 





- — a 4 a' 3 [ — ] A fj. (1-7 ^ 2 ), 
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2 





12 


U>| -J Co| 




B 4 = k 30 ^ a ' 4 (t) k " 3fJ) 

+ k 31 V ^ a>4 (f) ^ <*- 2 -3/i 2 ), 

B 5 = k 30 °- 3 a ' A ^7 j ^( 3 ^ 2 -M 2 ) 

+ k 31 - — a 4 a ' 4 ^ /x v (3 \ 2 - fj. 2 ) , 
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B e = +y k 3i V a2a<4 (?) (x2 "^ 2) 

+ k 30 i— a 3 a' 4 ^ ^ (\ 2 - ,u 2 ) 

- I k 31 a 4 a' 4 ( | y (X 2 - M 2 ) (1 ~ 7 v 2 ), 

„ 1 , a ' 4 2 <4 ( a \ 3 v 

B 7 = +y k 31 — a a [j]^ 

+ k 3Q — — a 3 a' 4 ^ kfiv 

(±YV(l-7v’), (S’) 

It can be checked easily that in outer space A . and B ) satisfy the Laplace 
differential equation. On the surface of the Earth, using 

£_ = 1 + e 2 v 2 + • • • • , 
r 2 

Aj and B. become, respectively, and bj . These expansions represent a 
generalization of the expansions obtained previously (Musen and Estes, 1971) 
where the assumption is made that k 2 and k 3 are constants. 

The terms of higher order in e 2 are not of any importance from the ob- 
servational standpoint and they are therefore omitted. 
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Substituting (5) and (5’) into (1) instead of a. i (j = 0, 1, .... 4) and bj (j = 
1, 2, . . . , 7) and making use of 

GM = n 2 a 3 


we obtain the following representation of the tidal disturbing function in outer space: 


“■S --(I- 

+ n2a2 ¥ aa ' 3 ^l (-fr) (7) V 1 


2 2 m 
+ rr a - 4 — 

M 


r * 2 “' 3 (?0 3 (f) 3v ’ 




(6) 


V 2 is the dominant part. It defines the so-called "main problem." It is of the 

• ■* . .s' , 

second degree in k , f± and v . We have 


v a = p;p„+ p;p, *p;p, + p; p 3 + p ; p 4 . 


( 7 ) 


where 







b ;- 
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P 1 “ ( k 20 “ J **22] ^ ~ Y4 ^ A 1 + J k 31 a ~ B 6> 


Pj = k 2 - fJ. 2 , 


P 2 _ ( k 20 "y k 22^ 14 A 2 + y k 31 a B 7* 


P, = X/i, 


P 3 


= ( k 2 o+y k ”) (^if 62 ) A 3 ~y k 3i a 'y7 B l ( 


p 3 =^V, 


p ; 


- ^ k 20 y k 22^ ( 1- l4 £2 ) A 4~y k 3i a 7? B 2 • 


P 4 =? M v • 


(V) 


The term V, is of the third degree in K, ^ and v . It contains the k 30 part 


and the influence of the deviation of k 2 from a constant. 


v 3 = s; s t + s' s 2 + s' s 3 + s; s 4 + s; s 5 + s; s 6 + s; s 7 , 


where 


(8) 


Sj — k 21 A 3 + k 30 a r# Bj , 


Sj = \ (1 - 5 v 2 ) , 
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C' - _ k 
& 2 ~ 5 K 21 


a: 


+ k 30 a ' 


B', 


S 2 = ii (1 - 5v 2 ), 


S 


3 



A' 


'■30 




S 3 = v (3 - 5 v 2 ) , 


s ; = k 3 o a 'A-B;, 


s 4 - k (X. 2 - 3/x 2 ), 


s 5 (3\ 2 - M 2 ) , 

s ;^ k 2 ! a; + k 30 a'^ B ;, 
s 6 = V (\ 2 - M 2 ) . 

S y = k 21 A 2 + k 3Q a — 7 B ? , 


S 7 = X./J.V. 


( 8 ») 
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V 4 contains the influence of the oblateness and also the influence of the 
deviation of k 2 and k 3 from constants. It is of the fourth degree in k, fj. and v. 
We have 


v 4 = t; Tl + t' t 2 + r 3 t 3 + t; t 4 + t' t s + t; t 6 + t; t 7 . 


O) 


where 



T 4 = (\ 2 - 3 /x 2 ) , 
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T g = fj. v (3 X? - /j. 2 ) , 



T 6 = (\ 2 -/z 2 ) (1-7 v 2 ) , 



T 7 =^(1-7v 2 ). (9') 

The factor (a'/r 1 ) 3 A' 0 attached to GM/r in the first term in (1) contains the 
purely periodic terms with the arguments i, V , F, D and r and a constant 
term. Consequently, the influence of the oblateness and the deviation of k 2 
from a constant on the tidal potential may affect the accurate determination of 
GM. 

V, = -±a^ + I(a;\ + a;m) 

The V 1 terms are not of any importance in the computation of the tidal 
effect in the motion of an artifical satellite because they disappear in the proces 
of elimination of the short period terms by averaging the disturbing function 
over the orbit of the satellite. 
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Substituting 


A. = + cos Ac os (f + o>) - sinftsin (f + co) cos i, 


/i = + sin flcos (f +<u) f cosAsin (f +w) cos i, 


v - sin (f +c<j) sini (10) 

into the expressions (7*) for P 0 , P 1 , . . . , we obtain 


P n = - — + — cos' 


3 

+ — sin 2 i cos (2f + 2 cS), 


= + -i- cos 2 A s in 2 i +^cos 2 A (1 + cos 2 i) cos (2f + 2 oS) 


- sin 2 Acos i s in (2 f + 2 co). 


+ isin2Asin 2 i + isin2A(l + cos 2 i)cos (2f +2a>) 


+ — cos 2 ftcos i sin ( 2 f + 2 w), 


P 3 — sin A sin i cos i + -^sinAsin i cos i cos (2f + 2oj ) 

2 2 


+ — cosAsin i sin (2f + 2 a), 


P 4 - + — cosAcos i sin i -icosftcos i sin i cos ( 2 f + 2 co) 


+ ySinAsin i s in (2 f + 2w) 
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and V 2 takes the form 

V 2 = V 20 + v 22 } cos ( 2 f + 2 w) + V<*> sin (2 f + 2 «). 

where 

v 2° = (--J + | cos2 ^ K ao ) - J K at ) sin 2i + ^ K 22 ) sin2 ( 11 ) 

V 22 } = +-| K 20 ) sin2 1 + ^ K 21 } sin 2i + J (! + cos 2 i) K^>, (12) 

V 22 } = + sin i - K 22 } COS l < (13) 


where we set: 

K< c ) - P' 

"■2 0 _ r 0 ’ 

^21 * = + Pj s in ft - cos SI, 

K 2 2 ^ = +PJ cos 2 ft + i P 2 sin 2 SI, 

K21 * - + P3 cos A + P4 s in fb, 

K 22 } = + P 1 sin 2ft- Ip' cos 2 ft • 

We also need the expansion of 

u n a n . . n 

n • V U = + — — sin 1 s in ft - — — sin 1 cos ft + - — cos 1, 

a A. o /j. a v 
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because it enters into the differential equations for the variation of SI and i . 
From (7*) and making use of (10) we deduce: 

R • V Q P 0 = - 3 sin 2i S i n (f + o>), 

R • V Q P 1 = + 2 sin 2ftsin i cos (f + qj) + cos 2 fi s in 2i sin (f + oS), 


R • V Q P 2 = - cos 2 ilsin i cos (f + oS) + -^ s in 2 ft s in 2i sin (f + w), 


R • V Q P 3 = + cos cos i cos (f +o>) - s in cos 2i sin (f + co), 
R • V Q P 4 = + sin Acos i cos (f + oS) + cos ft cos 2i sin (f +o>), 
and as a result: 


R * V Q V 2 = + V<‘> cos (f + oS) + V<0 sin (f + a>). 


where: 


V<«) = + cos i + 2K£| ) sin i, 


v( s ) 
21 


- sin 2 i - cos 2 i + K<*> sin 2 i. 


(14) 

(15) 


In a similar manner, substituting (10) into (8') we deduce: 


S 3 = - i(l-5 cos 2 i) cos Acos (f+w) + i(ll-15cos 2 i)cosisin<flsin(f+o)) 


+ ~sin 2 icosftcos (3f + 3w) - isin 2 icos i sinftsin(3f + 3w), 
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4*|00 4*1 4*1 Co 4*1 Cn 4*| Co 4*| Cn 4*|i-* 


(1 - 5 cos 2 i) s in A cos (f +cj) -- (11 - 15 cos 2 i) cos i cosftsin(f + w) 


sin 2 i sinftcos (3f + 3a >) + — s in 2 i cos i cosftsin (3f + 3 oj), 


(1 - 5 cos 2 i) s in i sin(f+oj) 


s in 3 i s in (3f + 3o>), 


sin 2 i cos 3 Ac os (f + oj) 


sin 2 i cos i sin3Asin(f+o>) 

4 


(1 + 3 cos 2 i)cos3Acos (3f + 3o>) - -i-(3+cos 2 


i)cos i sin3Asin(3f + 3co), 


s in 2 i s in 3 ftcos (f + o>) + — sin 2 


i cos i cos 3Asin (f + o>) 


+ i_Cl + 3cos 2 i)sin 3 Ac os (3f + 3 oj) + — (3 + cos 2 i) cos i cos 3 As i n(3 f + 3o>), 

4 4 


S 6 = - -isin2isin 2Acos(f + oj) + A-(l-3cos 2 i)sinicos 2As in (f + o>) 


+ i s in2i sin 2Acos (3 f + 3co) + i(l + cos 2 i)sinicos 2 As in(3f + 3w), 


S ? = + isin2icos 2Acos (f +w) + — (l-3cos 2 i) sin i sin 2Asin(f + co) 
8 8 


- — s in 2 i cos 2 Ac os (3 f + 3 oj) + -L(l + cos 2 i)sinisin2Asin(3f + 3 co). 

8 8 
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Substituting these values of S 1 , S 2 , . . . , S 7 into V 3 we obtain the following 
representation of the third part of the tidal'disturbing function 0 3 and of R • V Q fi 3 : 




a 3 a ' 3 



( 16 ) 


• {V^> cos (f +cd) + > sin (f +<o) + cos (3f + 3w) + sin (3f + 3w)}, 

R • V Q = n 2 a 2 • — • a 3 a' 3 \ ^ • (17) 

03 M \r / \ r / 

* ^ V 30 + V 32^ cos C 2 ^ + 2aj ) + v 32 5 sin ( 2f + 2w ))> 

where 

\ (1 - 5 cos 2 i) K(;> + 1 K< c 2 ) Sin 2 i + | K< c 3 > sin 2 i (18) 

^ = -2 (1 -5 cos 2 |i) sin i 

-1(11 -15 cos 2 i) K<;> cos i +1 (1 -3 cos 2 i) K<*> sin i, 

■5 

+ — sin 2 i cos i * (19) 

= +|-K^> sin 2 i -l^^ sin 2i + 1 (1 + 3 cos 2 i) K<*>, (20) 
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V<* > = + — K{ S J sin 3 i + - ) sin 2 icosi + ^(l+cos 2 i) 5 sin i 

o3 4 3U ^ 31 4 

i 

+ 1(3 + cos 2 i) K<“>, (21) 

V 3Q = - 1 cos i • (3 - 5 cos 2 i) 1 (1 - 5 cos 2 i) sin i 

+ | sin 2 i cos i • K<*> -|k< 5> sin 3 i, (22) 

^ sin 2 i cos i -l(l-3cos 2 i)K^ sin i 

2.2 

r * ■ •» 11 

-.1(1-3 cos 2 i) cos i --l(l+cos 2 i) sin i, (23) 

V<*> = -iK<;> sin 2 i + K<'> cos 2i +1 k<|> sin 2i, (24) 

2 2 

and 

= + S[ cos SI + Sj sin SI, 

K$5> = - s: sin 2SI + 1 Si cos 2A, 

32 6 2 7 

cos 3A + Sg sin 3A, (25) 
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k$;> = s;. 

Kjj > = - Sj sin A + Sj cos A. 

K<*> = + s; cos 2 A +^-S; sin 2A, 

^ 33 ^ = " ^4 S ^ n ^A + ^5 COS (26) 

In a similar manner we deduce: 

V 4 = V 40 + V<^ cos (2 f + 2oS) + V< s 2 > sin (2 f + 2w) 

+ ^ 44 ^ cos (4f + 4w) + V<* ^ sin (4 f + 4w), (27) 

where 

V 40 = + - (3-30 cos 2 i + 35 cos 4 i) K<'> 

8 

+ JL sin 2 i • (3 — 7 cos 2 i) - — sin 2 i (1 - 7 cos 2 i) 

16 8 

- — sin 3 i cos i • ^ , (28) 

8 


V 42 } = " f sin2 i (1 - 7 COS 2 i) 

— L sin 2 i • (4 - 7 cos 2 i) 

12 41 

+ .l(l_6cos 2 i+7 cos 4 i) ^ 


1 

cos 

2 


3 i sin i 


K (c) 

^43 ’ 


(29) 
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v( s ) - _ — (l-7cos 2 i)sin i • K< 8) 
v 42 12 ^ • 

+ 1(5-7 cos 2 i) cos i -K< 8) 

+ i_(l -3 cos 2 i) sin i • K^, (30) 


V<=> = + 11 sin 4 i *K<®> sin 3 i cos i • K<* } 

44 g 40 24 


where 


+ | (1 +COS 2 i)sin 2 i -K< c 2 > +^-(3+cos 2 i) -sin2i *K<'\ 


( 31 ) 


V ( s ) = + 2 - sin 3 i * K<!> - \ sin 2 i cos i *K< 8) 

44 94 41 4 


24 


+ 1.(1 + 3 cos 2 i) sin i • K*j 8) , 

8 


( 32 ) 


R • V 0 v 4 = V<' } cos (f +«) + V< 8) sin (f + w) 

+ V^ } cos ( 3 f + 3 w) + V< 8) sin ( 3 f + 3 w) , ( 33 ) 


v( c ) 

v 4 1 


.1.(3 -7 cos 2 i) cos i -K<»> -|sin i (1-7 cos 2 i) • K^ 8) 


Isin 2 i cos i *K< 8 >, 


(34) 
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^41 ^ = + 15 sin i cos i (3 - 7 cos 2 i) • 

- 1(3-27 cos 2 i + 28 cos 4 i) -K<J> 

- 3 sin i cos i (4 - 7 cos 2 i) • 

+ |-sin 2 i (1 -4 cos 2 i) • K<5>, (35) 

= + — sin 2 i cos i • ^ ^ sin i • (1 - 3 cos 2 i) • 

4 2 

- i (5 - 9 cos 2 i) cos i • ) , (36) 

= - 35 sin 3 i cos i • ^ sin 2 i (1-4 cos 2 i) • 

- 7 sin i cos 3 i • i (3 - 3 cos 2 i - 4 cos 4 i) * (37) 

and 

K(°) - + t' 

“40 _ + *3 * 

K<*> = + Tj sin ft - Tj cos ft, 

K 42 } = + T e cos 2ft + -|t' sin 2ft, 

• K<|> = + T; sin 3ft - T' s cos 3ft, (38) 
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^ 41 5 .= + Tj cos A + T' 2 sin A, 


K(s) 

^42 


+ T' sin 2 A Tl cos 2 A . 

‘6 2 7 


^ 43 ^ = + "*4 cos 3 A + Tj sin 3 A . 


(39) 


DIFFERENTIAL EQUATIONS FOR VARIATION OF CONSTANTS 

Only the long period tidal effects in the elements can be easily observed. 

To determine them, we make use of the classical equations for variation of ele- 
ments with the short period terms eliminated by averaging the right hand sides 
over the instantaneous orbit of the satellite: 


de 

dt 


na 2 e 


2 77 



(40) 



+ 2 sin 2 


1 dA 

2 dt" 


. dA 1 1 

sin 1 — = + ... ' — 

dt na 2 /TT^- 277 



sin (f +oj) dg, 


(42) 
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dL 2 1 f /l /? 2\ 

= - ■ r — — dg+(l-/l-e 2 ) 

dt na 2 277 J 0 Br 


dn , ^2 i d «< 


ri + 2/1 -e 2 sin 2 111! 
cTt 2 dt 


In our case the disturbing force can be conveniently represented in the form: 


vn = liir° + Vr° • V n O. 
or u 


V o n is the gradient of Cl with respect to r° and 


Vr° =1 (I - r° r°). 


From (44) and (45) we deduce: 


r 11 = r R • Vi! = R • V Q Cl 


Substituting (47) into (42) and (43) and taking 


, r 2 df 
dg = — • - 

a 2 /l - e 2 


into account, we deduce: 


= -— ■ — 

na 2 e 277 J n \ a / 3 f 
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1 


1 

2t7 


dv 

dt 


na 2 e 



B fl , 

- a — cos r + 
Br 



BO 

S 7 



df 


+ 2 sin 2 


1 d<n, , 

2 dt ’ 


( 49 ) 


. . dA 1 

sin l = + 

dt na 2 (1 - e 2 ) 


~ f UY (R-V 0 fi)sin(f +«)df, 

Jo 


(50) 


di _ 1 

dt na 2 (1 - e 2 ) 


277 



(R • V 0 0) cos (f +oj) df , 


(51) 


dL 

dt 


2 1 
• ■ 

na 2 vTTT 7 217 



df. 


+ (1 - Vl - e 2 ) — + 2 /l - e 2 sin 2 — 
v ' dt 2 dt 

In the main problem we have: 

Q = n 2 a 2 — a 2 a ' 3 ( lY [V 90 + V<| >cos (2f + 2co) + V<|>sin(2 f + 2a>) ] (53) 

2 M \t 7 \ r / . . 

and 

R • V Q = n 2 a 2 — a 2 a' 3 f f— ^ [V^ 'cos (f + a>) + ) sin(f + &)) ] (54) 

02 M \r' \r 21 21 
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where V 20 , V<'> , V<*> are given by the equations (11) - (13) and 
by (14) - (15). 

Substituting (53) and (54) into (48) - (52) we obtain, after easy transformations, 
the following expressions for the perturbations as caused byft 2 : 


dS e 

— — = 0 
dt 


d 8 2 7T 

dt 


3n a 2 a ' 3 
(1 -e 2 ) 2 



+ 2 sin 2 


^ d8 2 A 
2 dt 


sin i 


d§ 2 <a 

dt 


2 * 3 

n or a ° 


2(1 -e 2 ) 


2x2 


m • _ V< s > 

¥ Ir' 1 21 ’ 


d^i 

dt 


2(1 -e 2 ) 2 


m' 

¥ 



dS 2 L 

dt 


6n a 2 a' 3 
(1 - e 2 ) 3/2 


/ 

m 

¥ 



1 + V 1 - e 


dS 7r : dS 2 ft 

+ 2 /l - e 2 sin 2 

I dt 2 dt 


In limits of the "main problem" the eccentricity of the artificial satellite 
is not affected by the tidal perturbations. This result was obtained by Kozai 
(1965) using the differential equation for the general perturbations in satellite's 
eccentricity. 

For the perturbative effects caused by ft we obtain, making use of (16) - (17): 


d S e 

dt 


(1 -e 2 ) 2 



sin co - V^* 1 cos co) , 
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^a 77 na 3 a' 3 


3 


dt (1 - e 2 ) 3 e M 


— ■ (£.]■ • (1 +4e 2 ) (V<'> cos a> + V<*> sin co) 


j'dS.ft 

+ 2 sirr 4 - — , 

2 dt 


sin i 


dS 3^ na 


d t 


(i 


3a<3e . HL . (£\ 3 . (' 

_ e 2) 3 M \v') \ 


V, n sinw- - sina> + 1 ) cos a> 

30 2 2 


d ^3 * n a 3 a ' 3 e m' 


dt (i _ e 2 ) 3 M 


[T) 3 • ( V 3» 


cos co + 1 cos a> + 1 Vl® ^ sin a> ) , 

2 32 2 32 / 


d ^ 3 ^ 8na 3 a' 3 e m' (a'\ 3 .. r ( c \ ,.f s n \ 

= + • — • I — - I * (V^5 > cos co + V^® ' sin an 

dt (1 _ e 2 )5 /2 M \ r / 31 31 


dS,7T 


e “ ^3 " _ r. r . 2 i 3 

2 v 1 - e z sin . 


1 + /l - dt 


1 dS 3 A 

2 dt 


Si mil arly, for the effects caused by the part fi 4 of the disturbing function 
we have, making use of (27) and (33): 


d ^4 e 3na 4 a' 3 e m' ( a'\ 3 , ..( C ) • o \r(s) 0 \ 

+ ( + V 42 Sm 2w - V 42 COS 2 ")’ 


dt 


2(1- e 2 ) 3 M ^ r 


dS 4 7T _ na 

dt “ + 


2-^1 •^•(K ) 3 f + flO + — e 2 ^ 

_ e 2 \ 4 M V r / \ \ 2 / 


(i-o 


'40 


^ y + ^- e 2 j ( + V ^2 ^ cos 2 co + 5 sin 2 co) 


+ 2 sin 2 


1 dS 4 ^ 

2 dt ’ 
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dS 4 <ft 

sin i — = + 

dt 


na4a ' 3 # m' t /a' \ 3 / /. 

2(1- e 2 ) 4 M" | + V 


1 + pj v<°> 


+ i e 2 [ + (-V«;> + v£>) cos2ai + ( + V<‘> -V<§>) sin 2w] J 


d8. 


dt 


c) 


i = + "« 4 ° ,a -ill • (4) 3 l(l + le>) VC 

2(1- e 2 ) 4 M \ r / 2 ' 

+ -|e 2 [ + ( + v<c> + V(O) cos 2 u> + (V<;> + Vj->) Sin 2w]J. 


d§ 4 L _ , 10 n a 4 a' 3 . m| . (V\ 3 . 3^ y 

dt (i _ e 2 ) 7/ 2 M \r'/ (\ 2 / 40 


+ — e 2 cos 2w + ^ sin 2 oj) 


1 — 

) 1 + 


_ dS 7T • dS.ft 

e 4 _ n. 2 ■ 2 1 4 

— + 2 v 1 - e z sin 


/m 5 dt 


2 dt 


CONCLUSION 

The present theory of the tidal effects represents an improvement over the 
theory based on constancy of Love numbers. The differential equations for the 
tidal effects in the elements of an artificial satellite are given in the present 
article in a form convenient for numerical integration over a long interval of 
time. This form is also suitable for representing the tidal effects as trigonometric 
series with the arguments a>, A of the satellite and l , V , F, D, r of the Moon. 

The numerical integration can be performed using some simple quadrature 
formula, without resorting to a predictor-corrector system. The final goal is 
the determination of the Love numbers k 2 , k 3 , . . . which we assumed to be not 
constant, but latitude dependent. This assumption required the reformulation of 
the "main problem" of the theory of tidal effects. 
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In the present exposition, the main problem contains the influence of oblate- 
ness and the deviation of k 2 and kj from constants. The effects caused by k 2 
and k 3 (including their deviations from mean values) are the most significant 
ones. It is extremely doubtful that the effects caused by k , k , . . . are significant 
enough to be detected. 

The computation done in the previous work (Musen and Estes, 1971) shows 
that the tidal perturbations in satellite orbital elements can be a few seconds 
of arc. The deviation of k 2 and k 3 from constants will probably contribute 
effects on the order of O'.'l - l'.'O. 

Some guess concerning a future direction of work can be made. We might 

/ 

expect that with the accumulation of observational material the problem of separa- 
tion of the long period tidal effects caused by the solid Earth from the effects 
caused by oceanic tides, as well as the problem of the non-homogeneous structure 
of the Earth, will arise. 

COLLECTION OF FORMULAS 
Basic information: 

Given: 

1° the mean elements of the satellite 

a, oo Q ,si Q , e, i, and the mean mean motion n Q , 
and<n, 0 are referred to the mean equator and equinox of a given initial epoch 

t o 

2° Lunar coordinates 

x\ y\ z' 

referred to the same equator and equinox. 
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3° We have (Brouwer, 1959): 


\ 


w = % + n o w i (‘-V 
ft = + n 0 (t - t 0 ) 

where 

"i = + \ y ‘ 2 (- 1 + 5 ^ + § 2 7 ^ 2 [(_ 35 + 24 ^ + 25 ^ 

+ (90 - 192 v - 126 v 2 ) Q 2 + (385 + 360 v + 45 rf) 0 4 ] 


+ — y! [(21 -9v 2 ) + (- 270 + 126 v 2 ) Q 2 + (385 - 189 v 2 ) 0 4 ] , 

16 * 

<0,! = - 3 y' 2 9 + 1 y' 2 [(- 5 + 12 rj + 9 rf) 0 + (- 35 - 36 v - 5 rj 1 ) 1 9 2 ] 
+ |y; (5 - 3 7,2) (30 - 7 0 3 )] 


and 


7 ^ + 2 


1 j 2 r 2 


• Tj 


r A 


3 , J 4 R2 

8 .4 


n"8 


3 /fT" 


COS 1, 


4° 

r' = /x ' 2 + y ' 2 + z ' 2 

= x'/r', = y'/r', v' = z'/r' 

50 

D 0 ' = l-3i'' 2 (3-5/ 2 ) 

d '=\ ,2 - m ' 2 e ;=- V ( i - 5 / 2 ) 

1 1 r 
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DI=X'm' E' = — M' (1 - 5 V' 2 ) 


- _ a 


D' = X' V' E' = 


D: = m' v‘ E' = — X' u' v' 

4 4 r • ~ 


K--~ 

6 r 


^-v< (X' 2 - M ' 2 ) 


E' =^X* (X' 2 - 3 m' 2 ) 

a r 


V C3X' 2 -m' 2 ) 


Main problem: perturbations produced by fi 


a 20 = +|( k 20 + y k 2 2 ) f 1 


42 / 


^20 - + 28 h 3i a ' 


3 (,_ 2,_ \ I, 5 , 2 


a 2i = vf 2 ° "y" 22 / \ 1_ i4 e 


■) 


P 2 1 = + k il a ' 

v- 28 31 


a n/> = + 
22 x 


3 ( k 2 o' 7 k 2 2 ) ( 1 _ n e2 ) 


^22 = + y k 3i a ' 


a 23 = tt 24 = + 3 ( k 20 + 7 k 2 2 ) f 1 - P £2 


^23 ^24 


C 200 = - y + | COs2 1 C 21 1 = + COsi 


S 210 = “ 3 sin 2 1 


C 201 = Sin 2 1 


C 212 = + 2 S ^ n 1 S 211 = “ cos-2 i 


c = + 2. s in 2 i 
202 + 2 


S 212 = + Sin 2 i 
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•<I| OJ 



K = * “20 D o + ^ K 


K( c > = + p 
20 r o 


p; = + a 2i d; + /s 21 e- 


= + P3 sinfl-P^ cos A 


P 2 = + *22 D J + ^22 K 


K< 5 > = + p; cos 2 ft + Ipl sin 2 ft 
22 1 2 2 


P 3 = ♦ *23 P 3 + ^23 


= + P3 cos A + sin A 


p; = + * 24 d; + fi 2 , e' 


K<!> = + P.' sin 2 n -- 5 -PI cos 2 ft 
22 1 2 2 


^20 + *200 ^23 '* + *201 ^21* + *20 2 ^22 ' 


nr = + s 2io nr * s.i nr + s 21J nr 


v( c ) - 
V 21 _ 


+ C 2U ^21^ + C 212 ^22^ 


d 8 A 

s in i = + 

d t 


n a2 a ' 3 • — ■ (^.Y y(s) 

>2\2 M \V Vzl 


2(1- e 2 ) i 


d §2 1 = + n ^ a ' 3 ■ m'. /a'y* y(c) 

o2^2 M \ r '/ 21 


d t 


2 (1 - e 2 ) 2 


d t 


2 77 3 n a 2 a' 3 m' /a'\ 3 „ , i d 

= + ‘ — ' I — 1 V,„ + 2sm 2 - 


(1 _ e 2 ) 2 M Vr V 20 


2 d t 


d S 2 e 
d t 


= 0 


^ ^2 ^ 6 n a 2 a' 3 m' /a'\ 3 e d S 2 7T 

. = + •—•(—]• + 


d t 


(1 - e 2 ) 


2\3/2 M \r 


I '20 


i + /r ="? d 1 


, i d 8„ A 

+ 2 /l - e 2 sin 2 


2 d t 
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Perturbations produced by 0 3 : 


% = *5o k 2l 


^30 “ + 4 k 30 a 


a 31 “ a 32 - 5 k 21 


*31 = ^32 = + | k 30 a ' 


a„„ = + — k„, 
33 4 21 


^33 =■+ ^ k 30 


tt 34 = + 3 k 21 


^34 = + 1Sk 30 a ' 


^35 ^36 “ + g" k 30 a 


c 300 = - 2 ( 3 " 5 cos 1 ) cos 1 c 


3 20 


15 . 2 . 

= + — sin' 4 i cos i 


S 320 = 0 


3 2. 5 r 

C 30 1 = + T (1 ~ 5 cos i) sin i c 32J = (1-3 cos 2 i) sin i S 32 i = " — - s i n 2 i 

2 2 2 


C 30 2 = + | Sin2 1 COS *■ 


' 322 


- (1 - 3 cos 2 i) cos i 

2 


s„„„ = + cos 2 i 

322 


3 ■ 2 • 

C,„_ = - — sn-r 1 

303 2 


~323 


— (1 + cos 2 i) sin i 


S 323 ~ + 2 Sin 2 1 


C 310 = 0 


S 3 io = ~ ^ (1 ~ 5 cos 2 i) sin i 


C 311 = - 4 C 1 - 5 cos2 i) 


S 3 ii = _ ^ (11 - 15 cos 2 i) cos i 


C 312 = + ^ Sin 2 1 


S 312 = + — (1 - 3 cos 2 i) sin i 


3 . 2 . 

c,, = + — s in' 1 1 

313 4 


3 - 2 ■ 

s„,, = +— sin' 1 i cos i 
313 4 
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Fo- + a 30 D 0 + ^30 E 0 


K (c) = 0 

30 


f ; = + a 31 D' + P 3l E; K<‘> = + f; cos + FJ Sin a 


F 2 = + a 32 K + ^32 E 2 = - F 3 


K ( ,'> = -F; sin 2 a + If; cos 2 a 


F 3 = + a 33 D i + ^33 E 3 K iV = + F 5 C ° S 3A + F 6 Sin 3 A 


f; - + a 34 d 2 + /3 34 e 4 


= + F 0 


F 5 ~ 


+ ^35 E s 


= - f; sin A + F; cos A 




F: = +^ 36 E 6 = + F 3 cos 2A + -F; Sin 2 A 


Kl?) = - F' sin 3A + F' cos 3 A 


33 


V 30 “ + C 300 K 30 } + C 301 + C 302 K 32 } + C 303 K 33 } 

V 32 } = + C 320 K 30 } + C 321 K 31 } + C 322 QV + C 323 K 33 } 


v( s > = 

32 


+ S 321 E 3 1 ^ + S 322 E 32 ) + S 323 E 33^ 


v( c ) 

V 31 


+ C 3 1 1 E 3X^ + C 312 E 32^ + C 313 E 33^ 


V< 3 ;> = s 3 io *<,;> + s ..i K<;> + s 312 K<*> t s 3U «JJ> 


sin i 


d 8 3 A _ n a 3 a' 3 e _ m\ /a^\ 3 L s i n w _ Iv^ sin w + I V*®-* cos 

“TT-- + (1 _ e ^ M Vr ') \ 30 2 32 2 32 ) 


dS 3 i 


1 - | n a3 a e • m . • (v cos co + Iv< c J cos co + lv<’’ sinwl 

r “ „ " m W \ 30 2 2 / 


(1 - e 2 ) 


40 



d S 3 7 7 
d t 


•= + 


• — • 1 — ) (1+4 e 2 ) cos oj+ sin w) 

_ e 2)3 e M VrV v 31 31 


dS 3 e 

d t 


■ = + 


dS 3 L 

d t 


n a 3 a' 3 

„ t 

m 

(1 - e 2 ) 3 

e M 

2 sin 2 — 
2 

d S 3 A 

d t 

n a 3 a' 3 

. »' . (: 

(1 - e 2 ) 2 

M \i 

8 n a 3 a' 3 

e _ m' 

(1 _ e 2 ) 5/2 M 

e 

. dS 3 


31 


'31 


'31 


'31 


, : dU 

, + 2 /I - e 2 sin 2 — 

i + yrr^ d t 2 dt 


Perturbations produced by : 


a « + 56 \ 20 + 7 k ' 22 ) ^ " 280 k22 


/3 « = - 28 k ». 


45 /. 1 , \ 2 27 . 

a 41 “ a 42 " + J4 \ k 20 + f k 22j 14 


A,, = At, '= +-tf- k 3 l a 


22 *41 ' 42 


56 


a 


'43 


_ + 15 L 2 k U 2 - 9 k 

“ + 56 \ 20 7 k 22j 16 22 


/3._ = - II k„ a' 
43 28 31 


a. = +±£ fk_ ~?kj e 2 -ik 


15 


44 


14 


‘20 7 22 


22 


fi MA = - — k„ a' 
44 7 31 


a.. = a = 0 
45 46 


^45 = ^46 = + g k 3i a ' 
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C = + - (3 - 30 cos 2 i + 35 cos 4 i) 

400 g ' 

c = + 1_ (3 - 7 cos 2 i) sin 2 i 

401 + 16 v 


C 40 2 


-1(1 - 7 cos 2 

8 


i) sin 2 i 


3- 3 . 

c . = sir 1 cos 1 

403 g 


C 420 = - f^ 1 - 7 C ° S2 Sin2 1 S 420 = 0 

C 421 = * Y2 (4 " 7 COs2 ‘) Sin 2 1 S 421 = " f2 (1 ‘ 7 C ° S2 ^ Sin 1 

1 1 ? 

c 422 = + - (1 - 6 COS 2 i + 7 cos 4 i) s 422 = + — (5 - 7 cos' 1 ) cos 1 

2 ^ 


c 


423 


x 3 . 

_ cos° 1 s 1 n 1 

2 


423 


+ 1 ( 1 - 3 cos 2 
4 


i) sin i 


C 4 io = 0 S 410 ~ + ~2 (3 " 7 cos 2 i) sin 2 i 

c 4M = - 1(3 - 7 cos 2 i) cos i s 4U = - 1(3 - 27 cos 2 i + 28 cos 4 i) 

3 3 

C 412“ - ^ (I - ? cos 2 i)sini s 412 = -- (4-7cos 2 i)sin2i 

c = + — sin 2 i cos i s = + — (1 - 4 cos 2 i) sin 2 i 

413 4 413 4 v 


42 



0 


“430 


s.„„ ■= - 35 sin 3 i cos i 
430 


7 - 2 • 

c,„ = + — sin 4 l cos l 
431 4 


_ T _Z_ eir.2 

431 12 


+ — sin 2 i (1 - 


c 4 3 2 = +y(l - 3 cos 2 i) sin i s 


432 


= - 7 sin i cos 3 i 


c 433 = ~ ^ (5 - 9 COS 2 i) cos i s 433 =+i,(3-3cos 2 i 


G 0 a 40 E 0 + ^40 E 0 
G 1 = ' a 41 D 3 + ^41 E 1 


C 2 = *42 ^ 4 ^42 E 2 


G 3 = *43 K +^43 E 3 
= *44 D 2 + ^44 K 


G 5 ~ 


^45 E 5 


G 6 = A 6 E 6 


K (c) = G ' 

No u o 


=+Gj sin ft - G' cos A 
= + Gj cosft + Gj sin fi 


= + G' cos 2 SI + — G! sin 2 Tb 

42 3 2 4 

K<*> = + Gi sin2ft--G! cos 2ft 

42 3 2 4 


4 cos 2 i) 


4 cos 4 i) 


43 



= + G! sin 3 A - G' cos 3 A 

43 5 6 


= + G' cos 3 A + G ' sin 3 A 
43 S 6 


V = + c K (c) + c K( c) + c K (c) + c K( c) 

V 40 + 400 No + Noi *41 + 40 2 ^42 + 40 3 *43 


V( c ) = + c K^ C) + C K( c ) + C + c K( c ) 

V 42 + 420 *40 + M 21 *41 + 422 *42 + 423 *43 


y ( s ) - 

v 42 


+ s s ) + s + s K ( s ) 

+ 421 *41 + 422 *42 + 423 *43 


V( c ) = 

V 41 


+ c K( s ) + c + c K^ s) 

+ Sll *41 + L 412 *42 + 413 *43 


^41 + S 410 ^ 40 ^ + S 411 ^ 41 ^ + S 412 ^ 42 ^ + S 413 ^ 43 ^ 


y( s ) = 

V 43 


+ C 431 KV + C 432 K 42 } + C 433 K< 8 > 


^ 43 ^ + S 430 ^ 40 ^ + S 431 ^ 41 ^ + S 432 ^ 42 ^ + S 433 ^ 43 ^ 


d S A 


s in i 


V-= + ft ) 3 • {+ fl + — e 2 ) V<‘> 

t 2 ( 1 - e 2 ) 4 M U / l \ 2 / 41 


+ | e 2 [+ (- V 4 <"> + V 4 < 3 S >) cos 2 a) + (+ - V 4 < 3 C >) sin 2 w]} 


^ § 4 i n a 4 a * 3 /o'\3 f / o 

— J~=+JLU: • — •(*)■ \(i +±e 2 ) v< c > 

dt 2 ( 1 - e 2 ) 4 M \ r / l\ 2 I 41 


+ _ e 2 [+ (+ V 4 <'> + V< 3 C >) cos 2 c; + (+ + V 4 < 3 S >) sin 2 co] 


d 8 7 i 

4 

d t 



+ ^ + ^e 2 ) (+ V 42 ) cos 2 a> + V 4 ( 2 S > sin 2 co)| + 2 sin 2 -i ^ ... 
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d 8 e 

4 

d t 


3 n a 4 a' 3 _ m' 
2(1 -e 2 ) 3 M 



• (V 4 ( 2 C ) sin 2 co - ^ 


d s 4 l 

d t 


10 n a 4 a' 3 m' _ / a'\ 3 
(1 - e 2 ) 7/2 M \ r '/ 



1 + /I - e 2 ^ t 


d S 77 • d S 

+ 2 /l - e 2 sin 2 — 

2 d 


Indirect effects (Kozai, 1965): 


d S A 

= - (n fi. tan i) 8 i 

d t 


d Sj 7T 

— = (+ 5 sin i - tan i) n 


d S g L 

— = (+ 8 sin i - tan i) n, 

d t ' < 


cos 2 GO) 

e 2 ^ cos 2 co + V 4 < 2 *> sin2w) 

A • . • 


A x Si 
A 1 8i 
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